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We present a robust, ab initio method for addressing atom-light interactions and apply it to 
photoionization of argon. We use a close-coupling ansatz constructed on a multi-configurational 
Hartree-Fock description of localized states and B-spline expansions of the electron radial wave 
functions. In this implementation, the general many-electron problem can be tackled thanks to the 
use of the ATSP2K libraries [CPC 176 (2007) 559]. In the present contribution, we combine this 
method with exterior complex scaling, thereby allowing for the computation of the complex partial 
amplitudes that encode the whole dynamics of the photoionization process. The method is validated 
on the 3s3p 6 np series of resonances converging to the 3s extraction. Then, it is used for computing 
the energy dependent differential atomic delay between 3p and 3s photoemission, and agreement is 
found with the measurements of Guenot et al. [PRA 85 (2012) 053424]. The effect of the presence 
of resonances in the one-photon spectrum on photoionization delay measurements is studied. 

PACS numbers: 31.15.A-, 32.80.Aa, 32.80.Zb 



I. INTRODUCTION 

Since its discovery, the photoelectric effect has occu- 
pied a forefront position among the processes triggered 
by the interaction between matter and radiation. This 
is due to the relevance of photoionization for a range of 
naturally occurring phenomena, as for example those de- 
termining the opacity of astrophysical objects pQ, or lead- 
ing to radiation damage of biological systems [5] , as well 
as for many technological applications. A quantitative 
description of photoionization is required to understand 
and control these processes. Research in this direction 
has a long history [3J H] and generated a vast literature. 
The most recent development in this area includes the ad- 
vent of several new light sources providing shorter pulses, 
higher intensities or shorter wave lengths [5J |5J . During 
the last years, such sources have been utilized to gain 
knowledge also of the temporal aspects of photoioniza- 
tion 0|S]. 

A reliable description of the photoionization event re- 
quires a good representation of the atomic or molecular 
system during all the stages of its interaction with light. 
This includes the structure of the initial bound state, as 
well as of the final parent or target states, but also the 
coupling between the released photoelectron and the re- 
laxing parent ion. In this work, we present a tool that 
permits to treat many-body effects in bound and contin- 
uum states of many-electron atoms, irrespective of how 
the photoionization process is to be treated. Depend- 
ing on the intensity, pulse length and monochromatic- 
ity of the light the appropriate method may be found 
everywhere on the range from the standard one-photon 
time-independent formalism, to the explicit solution of 
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the time-dependent Schrodinger equation. Since the lat- 
ter approach quickly grows computationally heavy for 
all but the smallest systems, it is usually necessary to 
stick to a simplified treatment of correlation for larger 
systems. The procedure we outline here allows for a sys- 
tematic refinement of many-body effects and the unified 
approach for time-independent and time-dependent cal- 
culations permits further tests and verifications of the 
amount of correlation in the latter by comparison with 
the effect on the former. 

In the present study we validate the approach and ver- 
ify that the underlying approximations are justified by 
comparing photoabsorption cross sections for the argon 
atom, computed in the weak field limit, with existing ex- 
perimental and theoretical data. Furthermore, we show 
that the delay in the photoelectron ejection observed ex- 
perimentally [SI H] can be largely explained by the single 
photoionization scattering phases when corrected for the 
phase shift introduced in the experiment by the probe 
photon [TP] . 

For the description of the wave functions in the con- 
tinuum two different approaches dominate the litera- 
ture. One possibility is to match a numerical solu- 
tion, computed inside a inner spherical region, with 
the appropriate asymptotic solution in the outer region, 
whose analytical form is known. This is the basis of 
most scattering methods like, e.g., the R-matrix method, 
originally adapted for the photoionization problem by 
Burke and Taylor |llj and recently presented also in 
a time-dependent version [12], as well as the K-matrix 
method Ej ■ An alternative solution is to use com- 
plex scaling, where radial coordinates are rotated in the 
complex plane. In this case, the photoabsorption cross 
section can be directly computed as the imaginary part 
of the frequency-dependent polarizability j!5H17j . Com- 
plex scaling methods present some advantages as well as 
some disadvantages with respect to scattering methods. 
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On the one hand, with complex scaling, resonance states 
have a more transparent representation since they ap- 
pear directly as eigenstates of the complex-scaled Hamil- 
tonian with a negative imaginary part of the energy. On 
the other hand, the complex-scaled wave functions differ 
from the generalized eigenstates of the unsealed Hamil- 
tonian. Therefore, properties that can only be extracted 
from the ionized part of the wave function are not eas- 
ily obtainable. This problem is circumvented using ex- 
terior complex scaling |18j . where the radial coordinate 
is scaled only outside a sufficiently large region. In this 
method, the wave function in the inner unsealed region 
can be directly analysed. Exterior complex scaling has 
been extensively applied to electron and photon scatter- 
ing problems [TM2"2"] . see [23] for a review. 

The description of the target and parent-ion states is 
an atomic structure problem than can be treated with the 
tools targeting correlation in bound states, e.g. many- 
body perturbation theory, coupled-cluster, configuration 
interaction, etc., see e.g. [2"4Tf2"o] . Most of these ap- 
proaches can subsequently be combined with different 
representations of the continuum. Configuration inter- 
action, for example, can be used to either represent the 
atomic system in the inner R-matrix region or be em- 
ployed together with complex scaling. A method that 
is particularly convenient when aiming at problems in- 
volving single atomic continua, but non-trivial configura- 
tions, is the Multi- Configuration Hartree-Fock (MCHF) 
method [25] [57]. It is known to yield accurate corre- 
lation models for general bound states |25] while still 
limiting the number of configurations required for con- 
vergence. MCHF has often been used in combination 
with the R-matrix method, see e.g. Ref. [25], but also 
together with other approaches to describe the contin- 
uum (see e.g. Ref. [30!]). Here, we employ the close- 
coupling ansatz, originally introduced by Massey and 
Mohr [3TJ [35]; for the atomic wave function. The lo- 
calized part of the wave function is expressed in terms 
of MCHF configurations; our implementation is based 
on the MCHF atomic-structure package ATSP2K [33] . 
while the radial component that corresponds to the pho- 
toelectron is expanded on an extensive B-spline basis. Wc 
use exterior complex scaling for an explicit description of 
resonances and a good representation of the continuum. 

This paper is organized as follow. In Sec. [H] we outline 
the present implementation of the close-coupling ansatz, 



and in Sec. Ill the calculation of photoionization param- 
eters is discussed. To validate the method, we compare 
our results for the photoionization of the argon atom with 
other theoretical and experimental data in Sec. |IV[ Argon 
provides a good benchmark because it is comparatively 
light, so a non-relativistic description is appropriate, and 
its resonance profiles and angular distribution, which en- 
code the full information about the complex partial am- 
plitudes [34], have been the target of many experimental 
and theoretical studies [IH ESI I35H41"] . Fig. [I] is a dia- 
gram of the levels of neutral argon accessible through one- 
photon transition, showed on an energy scale correspond- 



FIG. 1. Diagram of the neutral argon levels and ionization 
thresholds accessible through one-photon transitions. Data 
are taken from the NIST atomic database [42]. The red level 
is the Ar ground state, and the blue states are singly ex- 
cited bound and autoionizing states. Black lines are ionization 
thresholds that are included in our analysis. Grey lines are 
other ionization thresholds which are neglected in the present 
paper (cf. Appendix \A\. 
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ing to photoionization from the ground state. Above 
the, 3s3p 6 threshold, we only mark the levels of Ar + 
in black and grey. In Sec. [iVj we focus our attention 
on the photon energy range between 26 eV and 30 eV 
which comprises the 3s3p 6 np autoionizing series, which 
entails strong configuration interaction. In Sec.[V]we use 
our approach to estimate the one-photon Wigner-like de- 
lay in the experimentally relevant region from 32 eV to 
42 eV. The importance of the presence of many reso- 
nances affecting the atomic delay is demonstrated. Fi- 



nally, in Sec. VI we draw our conclusions. 



II. THE CLOSE-COUPLING ANSATZ 

A. Exterior complex scaling 

In our implementation of the close-coupling ansatz, the 
states in the continuum are evaluated on an exterior com- 
plex scaled radial coordinate 



r 

Ro 



e l9 (r 



if r < Ro 
R ) i£r>R . 



(1) 



As mentioned earlier, this approach provides a good de- 
scription of the photo-absorption cross section in gen- 
eral, and of transiently-bound states in particular. In 
the present context, ECS has two further advantages. 



First, as detailed in Sec. II B we employ existing MCHF 
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codes [33 to describe the bound states. If i?o is large 
enough, the output of these programs can be directly 
used in conjunction with additional complex-scaled func- 
tions without any modification. In fact, this was the 
reason why Simon |18| . who wanted to treat molecules 
without having to scale the molecular potential, intro- 
duced ECS to start with. Second, with ECS we are able 
to analyze the ionized part of the wave function to deter- 
mine partial complex amplitudes. This will be discussed 
more in detail in Sec. Mil 

B. Multi-Configuration Hartree-Fock 

In the MCHF method, the atomic wave function is 
expanded on an orthonormal set of symmetry adapted 
linear combinations of Slater determinants [25j 127] . 
{$(7iL5)}, known as Configuration State Functions 
(CSF), 

* = J2cM^LS). (2) 

3 

The radial functions {P n i(r)}, which define the occupied 
and correlation spin-orbitals, and the {cj} mixing coef- 
ficients are then optimized to minimize the energy func- 
tional 

E[{P nl (r)},{cj}] = {*\H\V) (3) 

where H is the non-relativistic atomic Hamiltonian with 
infinite mass nucleus, under the orthogonality constraint 
for the radial orbitals {P n i(r)}: 

[Pni(r)P„>i(r)dr = 5 nn , . (4) 

We use the ATSP2K package [33] to solve the MCHF 
problem. 

The close-coupling expansion is defined by a set of N- 
electron target states ($ a ) coupled to an extensive set 
of one-particle states with defined orbital angular mo- 
mentum I. This group of (N + l)-particle states is fur- 
ther complemented by a set of localized states (x$)- The 
wave function of the atomic state before or after photo- 
absorption can be expressed as 

* = H c «X« + H Ccm ■A L s{*aPni(r N+1 )} (5) 

where a is the channel index, indicating target state a 
and photoelectron angular momentum /, and Als yields 
a symmetry adapted, antisymmetric (N + 1) - electron 
wave function. The iV-electron CSF are built from the 
orbitals obtained in the MCHF calculation, and we find 
the target states (\&a) as the linear combinations that di- 
agonalize the total ./V-electron Hamiltonian. We further 
use these orbitals to set up the (N + l)-electron CSF 
that form the Xi- All intermediate angular couplings 



compatible with the total symmetry under examination 
are included in the close-coupling expansion which, in 
this sense, is therefore complete. As in MCHF-based 
R-Matrix approaches [25], orbitals that are occupied in 
the target states, and which derive from a MCHF opti- 
mization, are kept frozen throughout the calculations. A 
discretized description of the continuum orbitals P n i(r) 
describing the photoelectron is assumed. The CI coeffi- 
cients, C{ and c an in Eq. ^ give the close- coupling ansatz 
a great flexibility. Even ifoccupied orbitals are optimized 
on the parent states, the second term of Eq. ^ is able to 
account for the orbital adjustment between the N- and 
N + 1-electron system. One drawback of this approach 
is that highly excited configurations or configurations in- 
volving correlation orbitals often appear as discrete states 
that are embedded in the continua, slightly above the 
physical states described by the model. These states ap- 
pear as artifacts of the calculation in the photoionization 
cross-section, so-called pseudo-resonances. 



C. Description of the photoelectron 

The photoelectron radial function basis {P n i}, cf. 
Eq. (pH), is expanded in B-splines [43] 

P*(rW) = 2>2J?(rW). (6) 

i 

The B-splines of order k are defined as piecewise polyno- 
mials of order k — 1 between pre-defined isolated points, 
called knots. In this work, we use k=7. With exterior 
complex scaling, it is convenient to let the knot sequence 
follow the scaled r in Eq. (fll), i.e. [44 . B-splines 
can be used to form effectively complete basis sets for 
the description of localized electron partial waves pack- 
ets, in particular in atomic problems [151 HB]. Since B- 
splines are non-orthogonal, however, they cannot be di- 
rectly used in combination with ATSP2K. The expansion 
coefficients in Eq. ^ are obtained through the following 
procedure. First the MCHF orbitals, computed on a ra- 
dial grid using the mchf program of ATSP2K, are fitted 
with B-splines using the w_bsw program of the B-spline 
R-matrix (BSR) package [57]. The exterior complex scal- 
ing radius, Rq, is chosen to be larger than the radius at 
which the MCHF orbitals numerically vanish, therefore, 
the MCHF orbitals are expressed in the region where r is 
unsealed. Second, for each orbital obtained at the MCHF 
stage, the B-spline corresponding to the largest mixing 
coefficient in P n i(r) is removed. Finally, the remaining 
B-splines are orthogonalized to the MCHF orbitals with 
the Gram-Schmidt algorithm. When combined, the oc- 
cupied P n i{r) and the orthogonalized B-spline set form 
an orthonormal basis equivalent to the original B-spline 
radial basis. The B-splines which do not vanish at Rq 
are complex [23] ; as a consequence, the final basis set is 
complex as well. 
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III. ATOMIC PHOTOIONIZATION USING 
EXTERIOR COMPLEX SCALING 

The correlated many- electron basis set that we use to 
describe the system before and after photo-absorption is 
obtained by diagonalizing the full complex scaled Hamil- 
tonian H e projected on the space spanned by the ansatz 
in Eq. The resulting eigenenergies E^, corresponding 
to the eigenvectors | k ), are in general complex. Those 
corresponding to bound states have vanishing imaginary 
part of the energy; some complex eigenvalues, which cor- 
respond to resonances, are largely independent of the 
rotation angle; finally, the remaining eigenvalues corre- 
spond to the unstructured continuum. 



A. The total cross-section 

The single-photon absorption cross-section can be ex- 
pressed in terms of the dynamic polarizability of the 
ground state |0) of the atom, 



4-7T U 
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Huj) (e • R) | 0) ,(7) 



where R = Y,j r j> G o( E ) = (E - H + iO+y 1 is the 
Green's function of the field-free hamiltonian, hut is the 
photon energy, e is the direction of the laser polarization. 

With complex scaling, the Green function is con- 
structed as a sum over the many-electron basis set com- 
ponents as [E3 [To] 
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where | k e ) and (k e \ are right and left eigenvectors of the 
complex scaled Hamiltonian H e , with complex eigenval- 
ues E k . For the field- free Hamiltonian H e , the left eigen- 
vectors are simply the transpose of the right eigenvectors, 
see e.g. @S]. The sum over | k e ) goes over all eigenstates 
to H s including both resonances and continuum states. 
We note that since resonances appear directly as eigen- 
states of H e the effect of a resonance on the photoion- 
ization process is accordingly accounted for by a single 
term in Eq. (l8|). 



B. The Fano profile 

As was shown by Fano [49] the cross section in the 
vicinity of an isolated resonance, situated at E r (rel- 
ative to Eq) and with autoionization width T, can be 



parametrized as 



o-ie 



(<?• 



1 



(9) 



where q is the so-called asymmetry parameter and e is 
defined as 



hio — E r 

r/2 



(10) 



while a a and aj, are slowly varying as a function of e. 

Looking at the contribution from a specific resonance 
in the sum on the right-hand side of Eq. ([8]), we can 
indeed put each such contribution in the form of Eq. (|9| . 
The derivation is detailed in Ref. |16| . here we just give 
the results. Each resonance, k. contributes with 
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where the reduced energy e k of the resonance |fe) is de- 
fined as in Eq. ( 10 ) , with E k — Re[J5fc] — E being the po- 
sition of the resonance, and = — 21m[Ek\. The asym- 
metry parameter q k , but also and cr^ in , are here con- 
stants and are determined by the resonance parameters 
and by the complex dipole matrix element between the 
ground state and the considered resonant state, the nom- 
inator in the term in brackets in Eq. (J8|. Labelling it 

Rk + Hk = (k e \e ■ R (e) |0) 2 [IS], the three constants are 
found to be 
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(12) 
(13) 
(14) 

(15) 



If the sum over all the other contributions in Eq. Q, from 
other resonances as well as from the smooth continuum, 
just give slowly varying contributions in the vicinity of 
the resonance k, then it is clear that the total cross sec- 
tion is well described by Eq. ([9|. This condition should 
be fulfilled if the resonance is far from other resonances 
and from thresholds. 



The resonance parameters, E r , T and q in Eq. (12) 



are often obtained from experiment through a fit to a 
Fano profile. For validation purposes, we compare these 



with our calculated values in Sec. IV Another parameter 



that is also obtainable from experiments is the coefficient 

p 2 m, 



(16) 
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Here a a can be calculated directly from Eq. ([13]), while 
(Ji, is not given by a m i n alone but from the full sum over 
states k in Eq. ([8]). However, in the limit that Eq. ^ is 
valid, <7f, can be obtained from Eq. (18]) at e = —q. 



C. Partial cross-sections 

Eq. Q gives a convenient way to obtain the total cross 
section. For the extraction of the partial cross sections we 
analyze the many-electron, one-photon perturbed wave 
function, H(u;, Tx, T2, ■ . ■ , rjy), which is a function of the 
photon energy as well as of the coordinates of all the 
particles in the system. The perturbed wave function 
solves the following inhomogeneous differential equation 



(H - E - hw)^(uj, rx, r 2 , . . . , rjv+i) 
-(e-R) * (ri,r 2 , 



,r N+1 ),(17) 



where the source term on the right-hand side describes 
the excitation from the ground state by the dipole op- 
erator acting on all electrons. The concept of perturbed 
wave functions was introduced by Sternheimer |51j . and 
it remains a useful tool also today, with novel applications 
emerging in the field of attosecond science [TO] • Here it is 
obtained through direct summation over our many-body 
basis as 



(k e 
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which represents the full many-electron wave packet 
formed by absorption of one photon from the initial state. 
To extract partial physical observables, S is projected 
on a specific channel, a, defined by a bound iV-electron 
state of the parent ion and a fixed angular momentum 
of the outgoing electron. As a result, we obtain a ra- 
dial perturbed wave function (PWF), denoted p s a {ui,r), 
which is a superposition of the radial functions P n i{r) 
introduced in Eq. ([s]) and Q . Fig. [2] shows the complex 
and imaginary parts of p^(uj,r). For r < R , p 6 a {u,r) is 
equal to the radial perturbed wave function without com- 
plex scaling, while in the complex scaled region, r > Rq, 
it is artificially damped. As a second step, the complex 
amplitude A a (cu) of the energy normalized PWF is ob- 
tained by fitting p a (u>,r) to a linear combination of the 
energy-normalized regular and irregular Gf Coulomb 
functions 

Re[p a (w,r)] = Im[Ji Q (a;)]F, B (-l/fc, Ax) 

-Re[A a (w)]Gf(-l/fe ) *r),(19) 

where, here, k is the photoelectron wave number and E 
is the corresponding energy. The fit is carried out for a 
radial distance r < R so that it is not affected by the 
complex scaling, but large enough that the asymptotic 
behavior is established. In practice, we compute 5 , and 
A a , for a specific M magnetic quantum number. In the 
following, however, observables extracted from the com- 
plex amplitudes are implicitly averaged over orientations 



except otherwise stated. The absolute square of the am- 
plitudes is related to the electron flux and, hence, to the 
partial cross section, as 
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Note that the amplitude (and the cross section) is that 
of a photoelectron in channel a. The spin-averaged, an- 
gular dependent photoelectron amplitude, with given M 
magnetic quantum number, is obtained after multiplying 
it by spherical harmonics coupled by Clebsch-Gordan co- 
efficcnts as 



Ac,m(u) Y lm(0, <p) (L a (M - m); lm\L a M), (21) 



where L a is the angular momentum of the A^-electron 
system and L a that of the total (N + 1) - system. 
Expression (21 ) can be used to obtain the angular depen- 



dence of photoelectrons ionized by polarized light 
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where P 2 is the second order Legendre polynomial and <f) 
is the angle between the polarization direction, and the 
direction of the photoelectron. The parameter /3(u)) is 
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where 0(a; (3) is a purely angular factor given by 

2(2/ Q + l)(2^ + l)(2i Q + l)(2L p + 1) 
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where the index refers to the initial state, and a — 
(a,l a ) and j3 — (b,lp) refer to the two interfering chan- 
nels. 



D. Extraction of resonance properties from the 
partial cross section 



Properties like p 2 , defined in Eq. ( 16 1, can be extracted 



from the partial channels as well as from the total cross 
section. For this we follow Kabachnik and Sazhina [M] 
and Sorensen et al. |36) . and express the amplitude of 
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FIG. 2. Radial perturbed function for the 3p — > ed channel, 
in length gauge, for e = — q in the vicinity of the 3s3p 6 4p 
resonance. The red and dashed blue curves are its real and 
imaginary part, respectively. 




the PWF as the combination of a resonant and a 
decoupled A® complex amplitude 



A a (e)=F(e)A$ + AZ(e) 

H £ ) = -v— = f-'TTl: 

i + e 1 + e z 



(25) 
(26) 



where e is defined as in Eq. (10). If the decay of the 



bound component of the resonance is treated to lowest 
order in the Coulomb interaction, these amplitudes have 
the same phase up to a tt shift [3H [SU] ■ We show in Sec- 
tion IV C that it remains valid in the case of the 3s3p 6 np 



resonances of Ar, and that the effect of the higher order 
contributions, present in our approach, is small. From 
A% and A% (e) , we can calculate p 2 [3U] as 
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where we have used the fact that the non-resonant back- 
ground cross section, o~b, is obtained from the sum of the 
squares of the decoupled amplitudes, A®. For the extrac- 
tion of A^ and A^ (e) , we suppose first that A^ can be 
considered a constant. Provided A% (e) is slowly varying, 
we are able to describe it around a certain e as a n-term 
Taylor expansion 



Then A?; can be obtained from n th 



derivative of Eq. ( 25 ) as 



Ai 



d n A a {e) (d n F{e 



de r ' 



dc h 



(28) 



The number of terms, n, used is decided from the con- 
vergence of A^; for a given e, and the stability of A^ 
when e is varied can be used for assessing the reliability 
of the procedure. Making a comparison between the val- 
ues obtained for p 2 from the total cross section Eq. ( fl6| , 
and from the partial amplitudes, the right-hand side of 
Eq. (27), is a test of the extraction of the partial cross 



TABLE I. List of mixing coefficients in the MCHF expansion 
of the lowest 2 S (3s3p 6 ) state of Ar+. 



Configuration state function (CSF) 



mixing coefficient 



3s3p° 

3s 2 3p 4 ( 1 D)3d 

3s3p 4 ( 1 S)3d 2 ( 1 S) 

3s3p 4 ( 3 P)[ 4 P]3d 2 ( 3 P) 

3s3p 4 ( 1 D)3d 2 ( 1 D) 

3s3p 5 [ 1 P]4p 

3s 2 3p 4 ( 1 S')4s 

3s3p 4 ( 3 P)[ 2 P] 3d 2 ( 3 P) 

3p 6 4s 

3s3p 5 [ 3 P]Ap 



0.83498723 
-0.52070520 

0.09817089 
-0.09655511 
-0.08964443 
-0.05536407 
-0.03427725 
-0.01888098 

0.00827551 
-0.00088679 



sections as well as of the Fano parametrization itself. 



IV. APPLICATION TO ARGON 3s -> np 
RESONANCE REGION 

A. Details of the calculations 

In the present section, we briefly outline the con- 
struction of the MCHF model. Our approach is sim- 
ilar to those followed by Berrah et al. [53] and by 
Burke and Taylor [TT]; the reader can find further de- 
tails in the original references. The spectroscopic or- 
bitals {Is, 2s, 2p, 3s, 3p} are optimized by performing a 
Hartree-Fock (HF) calculation on the 3p 5 2 P° state of 
Ar + . Then, the {3d, 4s,4p} correlation orbitals are opti- 
mized in a MCHF calculation on the lowest 2 S' state of 
Ar + including single excitations of the 3s3p 6 configura- 

tion, as well as the 3s3p 4 3rf complex. As illustrated in 
Table[TJ this state is characterized by strong configuration 
interaction. 

The set of CSF included in the diagonalization of the 
TV-electron Hamiltonian defining the 2 S parent states en- 
tering the close-coupling expansion (see Eq. is the 
same as the one used in the MCHF calculation. For 
the 2 P° parents, we include the single excitations of the 
3s 2 3p 5 , and the doubly excited CSF arising from the 

2 2 

3s 2 3p 3 3c? , 3p 5 3d and 3p 6 4p configurations. 

This model yields a threshold energy for the 3s pho- 
toionization of Argon of 29.73 eV, compared to the exper- 
imental value of 29.24 eV [S5]. The difference of 0.49 eV 
between theory and experiment is due to the lack of 
balance between the description of the ion and of the 
ground state of the neutral atom. Since this discrepancy 
is mostly reflected in a collective shift of all spectral fea- 
tures, we simply present all our results on a photon en- 
ergy scale shifted downward by 0.49 eV to compare with 
the experiment, rather than using any semi-empirical cor- 
rections. 

The set of localized configuration state functions X£ m 
Eq. ([5| comprises all single and double excitations from 
3s 2 3jr, all triple excitations with at least two 3d elec- 

trons, and the 3s 2 3p 2 3<i Ap configuration, i.e., all those 
listed in Table 2 of Ref. \M. 



7 



FIG. 3. Top panel: Details of the complex energies of the 
field-free hamiltonian of argon calculated for the 1 P° sym- 
metry with 6 — 0.055 rad (red circles) and 0.11 rad (black 
squares). Bottom panel: Comparison of the total calculated 
cross sections (thick lines), in length (full red curve) and ve- 
locity (dashed blue curve) gauges. The black dots are the 
experimental values of Ref. |29| . in arbitrary units, scaled as 
in the original paper. The thin lines are the cross sections 
computed by including only the 3s3p 6 np resonances in the 
basis set. 
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The final close-coupling expansion using B-splines is 
obtained as explained in Section [Il B|IIC| For small radii, 
a dense B-spline basis is necessary to describe correctly 
the inner orbitals. A 0.1/ Z an step size is used from 
the origin to 1/Z an. Then, the knot spacing increases 
exponentially up to 0.9 an. The knot sequence is then 
linear up to 80 an. It is generated by the W_BSW program 
of the BSR package [47]. To represent high-n Rydberg 
states, we extend the node set with an exponential se- 
quence of 100 nodes from 80 to R = 400 an. The final 
basis contains 249 B-splines. We ascertain that the re- 
sults are stable with respect to changes of the complex 
scaling parameters by varying R , the radial distance at 
which exterior complex scaling begins, from 50 to 80 an, 
and by doubling the complex scaling angle, 6 in Eq.|l]). 



B. Total photoionization cross-section and Fano 
profile parameters 

The upper panel of Fig. [3] shows part of the complex 
rotated energy spectrum for the l P° manifold, obtained 
from the diagonalization of the field-free Hamiltonian H e 
using two different complex rotation angles. As expected, 
the complex resonance energies are virtually independent 
of the rotation angle. For example, all the parameters 
of the 3s3p 6 np resonances agree to one ppm up to the 
n = 13 state. The pseudo-continuum energies, instead, 
follow a 6>-dependent path in the complex plane 48J . Also 
Rydberg states with n > 13, that do not fit in the com- 
putational box, show a f5-dependence. 



FIG. 4. Top panel: partial cross-sections for 3p ionization and 
total cross-section, in length gauge. Bottom panel: relative 
error of the sum of the above partial cross-sections and the 
total cross-section computed with Eq. (JsJ , for Ro = 65 an and 
O = 0.055 rad. 
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The lower panel of Fig. [3] shows total cross-sections 
obtained in length and velocity gauge. The thick lines 
correspond to the cross-sections obtained by including 
the complete sum over the basis set in Eq. ^ , while the 
thin lines are computed including only the 3s3p 6 np res- 
onance series by summation over the corresponding 
given in Eq. (111. The latter curves reproduce well the 



resonant structure of the cross-section, but the remaining 
non-resonant contributions are required to obtain the ab- 
solute value of the cross-section. The good agreement of 
the calculated total cross-section in the two gauges in the 
transmission window suggests that Ob in Eq. ([9]) is well 
described and that the discrepancy comes mostly from 

Table |n] lists the parameters for some of the 3s3p 6 np 
autoionizing states computed with our method, in length 
and velocity gauges, and compares them to Ref. [29],I38|. 
Our results are generally consistent with those obtained 
with the R-matrix method [29 with similar choices for 
the computational model and basis. The widths are 
slightly larger than the experimental ones and those ob- 
tained from the fits to the R-matrix calculation. The 
agreement of the q-parameter from the R-matrix calcu- 
lation is good, but both calculations show rather large 
deviations from the experiment. Ref. |29j attributes this 
discrepancy to the difficulty to accurately represent the 
non-resonant background in the close-coupling expan- 
sion. The marked gauge dependence of <r a discussed in 
the previous paragraph supports the fact than the con- 
tribution from the resonant channel is not yet converged. 



C. Partial cross sections and angular distribution 



As a test of the robustness of the method outlined in 



Sec. IIIC we compare the total cross-section computed 



TABLE II. Comparison of the resonance parameters calculated in this work with the ones fitted to R-matrix and experimental 
cross-sections for the 3s3p 6 np resonances. For each resonance, parameters computed in length gauge are given on the first line, 
and the ones computed in velocity gauge on the second line. 
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with Eq. Q to the sum over the partial amplitudes of the 
2>p photoemission channels computed with Eq. ( 20 ) . This 
is shown for the length gauge in Fig. [4] As can be seen 
in the lower panel, the relative error between the two ap- 
proaches is of the order of 10 ppm at most energies, which 
is close to our computational accuracy. For the high n 
resonances, spurious long range effects affect the channel 
perturbed radial functions. Therefore, their matching 
to Coulomb waves, see Eq. (19), is less accurate at the 



considered Rq, but the error remains consistently of the 
order of 0.1%. Above the 3s ionization threshold, the 
3s3p e ep channel accounts for the 10% discrepancy that 
can be seen in the for photon energies above 29.24 eV. 

We show in Fig. [4] that the window-shape of the res- 
onances comes from the dominating 3p — > ed channel, 
while the 3p — > es cross-section is instead enhanced on 
resonance. The partial amplitudes cannot generally be 
described by the same set of resonance shape parame- 
ters. Here a s has, for example, its maximum value at 
e = —q w 0.4, where the total cross section has a mini- 
mum. This is due to a destructive interference between 
Af and Af that suppresses the partial background am- 
plitude. This is illustrated in Fig. [5] where the markers 
correspond to the norm and argument of the calculated 
partial complex amplitudes in the vicinity of the 3s3p 6 4p 
resonance, while the curves are obtained from the sim- 
plified expression 



-1^1 



where rj 



1 for a — d, and rj = 



(29) 

-1 for a = s. 



Here, we drop the 3p 5 label in a for conciseness. The 
good agreement between the curves and data points il- 
lustrates the reliability of the extraction of A® (e) and 
A^ and that A® (e) is a slowly varying function. More 
interestingly, it means that A®(e) and A^ are in phase 
(Arg[A%(e)/A%} w 0) while Af (e) and Af are in phase 
opposition (Arg[Af (e)/Af ] w tt) over a wide range of e. 
As mentioned in Section [HI D[ these relations are rigor- 
ously valid if the decay of the bound component of the 
resonance is treated to lowest order in the Coulomb inter- 
action [Ml [SU]. However, while they hold with increased 
accuracy as the width of the resonance decreases, they 
are not exactly realized within more general models like 
ours. 

Another test of our procedure for extracting the A® (e) 
and A^ values, and in particular of A a (e), is the calcula- 
tion of p 2 . The values from the right-hand side of Eq. ( 27 ) 



agree with those obtained from the total cross section as 
explained in Sec.[ll]within a couple of parts per thousand. 
We compare the values obtained from Eq . ([27] ) to the fit- 
ted ones from Berrah et al. [25] in Table [TlpThe present 
results in length gauge agree with experiment up to the 
3s7p state, while the velocity gauge underestimates p 2 
slightly. The R-matrix fitted values are systematically 
below ours and the experiment, except in the case of the 
3s3p 6 4p resonance. 

Finally, Fig. [6] compares the calculated /3(w) for the 
3snp resonances, n — 4 to 9, to the experimental data 
of Berrah et al. [29]. Overall, the agreement is excel- 
lent. No convolution is applied on the theoretical results 
so that discrepancies appear close to the high n reso- 
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FIG. 5. Norm, in atomic units, and argument of the par- 
tial photoionization complex amplitudes versus the reduced 
energy of the 3s3p 6 4p resonance, (a) 3p — ¥ ed channel, (b) 
3p — > es channel. The crosses show 1 on 7 computed values. 
The lines are obtained from Eq. ( |29[ ). Red and blue indicates 
length or velocity gauge, respectively. 
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nances. The j3 parameter is, as expected, about 1.5 where 
the d channel dominates, far from the resonances (see 
Fig. |IJ), while the photoelectron angular distribution is 
anisotropic (f3 — 0) when the s channel dominates, i.e. in 
the vicinity of the minimum of each Fano structure. The 
structure in between those two extremes is well described 
in the calculations, which means that, despite the differ- 
ence between theory and experiment regarding g-values, 
the relative amplitude and phase between the s and d 
channels, leading to angular dependent interferences, is 
well represented. 



V. APPLICATION TO PHOTOIONIZATION 
DELAY MEASUREMENTS 

Photoionization by attosecond XUV pulses assisted by 
phase-locked IR laser-probe fields is a widely used ex- 
perimental scheme to gain temporal information about 



FIG. 6. Calculated beta parameter (/3) in length (full red line) 
and velocity (broken blue line) compared to experiment [29] 
(black dots). The vertical dashed line indicates a change of 
photon energy scale. 
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sub femtosecond pulse structures and dynamics [53] . This 
laser-assisted photoionization process is often interpreted 
as an attosecond streak-camera method, where the pho- 
toelectron will gain momentum from the instantaneous 
vector potential of the probe field depending on its ex- 
act time of ionization by the attosecond pulse. Assum- 
ing that the intensity of both pulses is weak, the cor- 
responding light-matter interaction can be treated us- 
ing lowest-order perturbation theory and the on-set of 
this streaking mechanism can be described with high fi- 
delity using only a limited number of complex transition 
amplitudes [54] . The delay-dependent probability mod- 
ulations of the photoelectron distributions are then di- 
rectly related to the phase differences of the two-photon, 
two-color matrix elements (XUV±IR). These so-called 
atomic phases in laser-assisted photoionization |55| . have 
been implemented for reliable characterization schemes 
of attosecond pulses [SS]. Recently, similar experimen- 
tal techniques have been extended to study the delay 
in photoionization from different initial states in neon 
[7] and argon atoms [HI H] ■ Theoretical work has shown 
that the total atomic delay, ta, measured in these ex- 
periments can be approximated as the sum of the one- 
photon Wigner-like delay in photoionization, tw, phis 
the so-called continuum-continuum delay, r cc , which is 
induced by the laser-probe field and the long-range ionic 
potential, so that [gl lMl 1ST] 

TA = T W + T cc . (30) 

In order to evaluate tw, we first apply the method pre- 
sented in Sec. |III C| in order to extract the partial wave 
amplitudes, A a (u) which contain the asymptotic phase- 
shifts, S a (oj) = Arg[A a (oj)], relative to the Coulombic 
phase-shifts, ct;(cj), where I is the angular momentum 
of the photoelectron in channel a. Then, using the to- 
tal asymptotic phase shift, rj a {uj) = er;(w) + 8 a {uj), we 
construct the Wigner-like delay for channel a as a finite- 



10 



difference derivative 



.(a) _ Va(u>) - r) a {u<) 



where the difference between the XUV frequencies that 
will interfere due to the influence of the probe field are 
separated by two IR probe photons, w> — uj < = 2lojr. 

The absorption of a laser-probe photon, required to 
make the transition from the intermediate momentum ki 
to the final momentum kf, induces a phase-shift denoted 
4> cc . Starting from the formalism presented in Ref. |10j . 

we use the soft-photon approximation, k^/2 3> loir, and the circles and squares, respectively, 
find the following expression for the phase-shift, using 
atomic units (H = 1, e = 1, m e = 1, 47reo = 1), 



FIG. 7. (a) Wigner-like delays for photoionization from 
argon. There is no perceptible gauge dependence for pho- 

(31) 

toionization from the 3p orbital. For the 3s photoionization, 
the results from velocity (blue broken) and length gauge (red 
dashed) are shown, (b) Wigner-like delay difference between 
photoionization from the 3p and 3s orbitals in argon. The thin 
curves are obtained with the full model that exhibit a pseudo- 
resonance at 39.4 eV. The thick curves are obtained with the 
truncated model for velocity (blue dashed) and length (red 
full) gauge, respectively. The experimental values of Ref. [5j 
are given with error bars. The experimental data point ob- 
tained using the r cc from Eq. ( 33 1 and Ref. [5S] are given by 



(j) cc (kf,UJ IR ) = arg 



i + 



k 3 

'(32) 

where Z = 1 is the charge of the target ion. Eq. (32) 
is accurate for photoelectrons with high kinetic energy. 
The continuum-continuum delay can then be written as 



z(k f ,aj IR ) 



4>cc(kf,UlR) 
OJIR 



(33) 



Alternatively this delay could be estimated with other 
methods such as the eikonal Volkov approximation [55] 
or classical-trajectory Monte-Carlo simulations [55] . 

If there is more than one channel contributing to a 
given final state, the two-photon amplitudes have to be 
added in a consistent way depending on the experimen- 
tal observable (TO] , Such a detailed analysis is beyond 
the scope of this paper, where we will show that the 
correlated Wigner-like delay Eq. (31 1 from the dominant 



channels is enough to explain the observations within the 
experimental uncertainty. 

Kliinder et al. [5J and Guenot et al. [9] have mea- 
sured the 3s — 3p atomic delay difference at XUV pho- 
ton energies of 34.1, 37.2 and 40.3 eV corresponding to 
the high-order harmonic sidebands 22, 24 and 26 for a 
generating field of 1.55 eV. The contribution from the 
continuum-continuum delay was then subtracted accord- 
ing to Eq. (30 1, in order to isolate the Wigner-like delay 



in photoionization. 

The first pseudo-resonance in our model lies at an en- 
ergy of 39.4 eV, which is inside the range of the exper- 
imental measurement, leading to an unphysical feature 
at 37.8 eV. In order to assess the reliability of the com- 
puted delays, we use a reduced close-coupling expansion 
by excluding the N- and (N+l)-electron CSF involving 
the 4s correlation orbital, which is responsible for this 
artifact. In this truncated model, the energy scale is 
shifted downward by 0.45 eV to match the Ar 3s photoe- 
mission threshold. Another pseudo-resonance remains at 
41.3 eV, affecting the Wigner-like delay at energies close 
the sideband 26, preventing comparison with the last ex- 
perimental value. 
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In Fig. [7] (a) we present the Wigner-like delay obtained 
with the reduced model. It is computed by Eq. (31 1 for 



channels leading to photoemission from the 3p and 3s 
orbitals. In Fig. [7] (b) we show the delay difference be- 
tween the two dominant channels, 3s3p 6 ep and 3p 5 ed, 
using the full and reduced close-coupling models. An 
exact determination of r cc is difficult and we, therefore, 
compare these curves to the original experimental data 
points for the Wigner-like delay of Ref. [9] together with 
the ones estimated with the alternative r cc from Eq. ( 33 ) 
and Ref. 



|58| . There is good agreement between our the- 
oretical model and the experimental points of sidebands 
22 and 24. 

Going beyond the present analysis requires to include 
the many shake-up thresholds opening at these ener- 
gies |52j , as shown in Fig. [I] As shown in experimental 
data [50] and in R- matrix calculations [5T], the 3s pho- 
toionization partial and the total cross-sections are highly 
structured by resonances associated to those thresholds. 
This suggests that the corresponding phases may be 
strongly altered as well. However, the widths of these res- 
onances are in general below 0.1 eV so that we can expect 
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that they do not drastically affect the measured delays 
obtained from an average over the sidebands. Still, as 
we will show below, some structure should be observable 
provided that the IR probe field is sufficiently narrow. 

The sideband amplitude can be obtained as the con- 
volution of the IR spectral envelope with the dipole el- 
ement between the perturbed wave function and the fi- 
nal state [55]. If, as before, we limit our discussion to 
Eq. (30), we can estimate the effect of the probe field 
bandwidth Tjh by extracting the Wigner-like delay from 
the complex amplitudes convoluted with a gaussian func- 
tion 



/ 41n2 



dx A a (x)exp 



41n2(w - 



r 2 

1 IR 



(34) 

In Fig. [SJ we present Wigner-like delay differences ob- 
tained with a model that includes a range of physical 
resonances associated to the thresholds in black in Fig. [T] 
(details are found in Appendix A| , after smoothing by 
the bandwidth of the probe field using Eq. (34). The 



top panel shows the results for a gaussian IR envelope 
with FWHM of 60 meV, corresponding to an IR field 
duration of 30 fs, typical of experiments using trains of 
attosecond pulses. The bottom panel corresponds to a 
FWHM of 0.6 eV, corresponding to a single IR cycle 
probe pulse, typical of experiments using single attosec- 
ond pulses. We show that the presence of resonances 
does not change radically the baseline of the relative de- 
lays, but rather superimpose a structure that preserves 
the previous theory-experiment agreement. Including the 
physical resonances in the model also pushes pseudo- 
resonances higher in the spectrum so that a comparison 
can be made with the experimental delay extracted at the 
sideband 26. Here, some limited agreement can still be 
found in the top panel, within the theoretical and exper- 
imental uncertainty. The Cooper minimum emphasized 
in Ref. [S], experimentally located at about 43 eV, is also 
located at higher photon energies in our calculations so 
that it does not affect the computed delays. We should 
also note that, above 39 eV, many more resonances could 
affect Tyy, corresponding to shake-up electrons in higher 
orbitals, even though they seem to be less important as 
far as the cross-section is concerned [60]. 

An important implication of these results is that a rich 
structure, originating from resonances, could be observed 
experimentally by scanning the energy of the high-order 
harmonics, see Fig. [8] (a). In order to observe such fea- 
tures, the bandwidth of the IR field should be sufficiently 
narrow. Otherwise, if a too short probe pulse is used, 
which is often the case in attosecond streaking experi- 
ments, these sharp energy-dependent features cannot be 
resolved. In this case, the resonances may affect the mea- 
sured delay smoothly, as is shown in Fig. [8] (b), where 
they cause an additional delay of up to several tens of at- 
toseconds as compared to the model without resonances 
in this energy region. 



FIG. 8. Relative Wigner-like delays using a model including 
some doubly excited resonances (see Appendix: A), (a) The 
probe field bandwidth is 60 meV and 0.6 eV, respectively, 
for the top and bottom panels. Length (continuous red) and 
velocity (dashed blue) gauge results are showed. The thicker, 
light curves correspond to the reduced model used in Fig. [7] 
Experimental points are showed with the same notation as in 
Fig.0 
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Such a streaking experiment has been performed 
in neon to extract a relative photoionization delay of 
21(5) as between the 2p and 2s electron extraction [7]. 
No theoretical analysis has been able reproduce this re- 
sult so far. Two of them include the effects of a probe field 
within an ab initio many-body framework. Full-fledged 
simulations of the measurement have been performed by 
Moore et al [63] , using time-dependent R-matrix theory, 
resulting in a relative delay of 10.2(1.3) as. Later, we have 
reproduced this result using two-photon matrix elements 
computed using many-body perturbation theory [64] ■ In 
that study, the validity of Eq. (30) is demonstrated. In 



both theoretical approaches the effects of transient dou- 
ble excitations are neglected. As in the case of Argon, 
the Ne + level structure at the relevant energy, i.e. a XUV 
energy of about 105 eV, is extremely rich [55]. Here we 
have seen that such effects might play an important role. 
In Neon, an extra complication is that the experimental 
photon energy is above several double ionization thresh- 
olds, and that triply excited states may also be impor- 
tant, as pointed out in Ref. [63]. This implies that it will 
be important to include a complete description of at least 
two active electrons, and possibly excited configurations 
of the grand-parent system. 
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VI. CONCLUSIONS AND OUTLOOK 

We have presented a practical method, based on the 
ATomic Structure Package (ATSP2K) [33], to obtain a 
general, many-electron atomic wave function in a close- 
coupling expansion. The structure of the localized part 
of the wave function is represented using the Multi- 
Configuration Hartree-Fock (MCHF) method, which pro- 
vides a versatile and flexible way to represent complex 
many-body effects with a limited number of configura- 
tions. By using B-splines to describe the radial compo- 
nents of the photoelectron in the allowed partial-wave 
channels, we obtain a multi-channel method that, for the 
description of one active electron in a box, is complete. 
This paper develops this approach in combination with 
the exterior complex scaling method of use, in particular, 
in the description of monochromatic wave packets. In ad- 
dition to evaluating the total cross-section with the stan- 
dard tools of complex scaling, we have extracted complex 
partial amplitudes that encodes the full information on 
the studied process. The consistency between the two 
approaches has been verified in detail. 

The method is validated by application to the pho- 
toionization of Argon close to the 3s ionization threshold. 
In this energy region the final state is affected by strong 
configuration interaction. We found good agreement of 
the line profiles of the 3s3p 6 np autoionizing states with 
previous R-matrix calculations and experiments. Also 
the computed photoelectron angular distribution is in ex- 
cellent agreement with the measurements. 

It has earlier been suggested that laser-assisted pho- 
toionization by attosecond pulses gives access, in a uni- 
versal way, to relative one-photon Wigner-like delays (see 
e.g. Refs. [ZHS]). So far, this interpretation of the obser- 
vations has remained uncertain since no ab initio cal- 
culation have been able to reproduce the experimental 
data, even though some encouraging results were ob- 
tained in the case of Argon [9] [64]. In this paper, we 
present a study that contributes to validate these claims, 
in the context of a weak probe field, by including all 
relevant correlation effects, and by showing that, in so 
doing, satisfactory agreement is obtained with the ob- 
served atomic delay shifted by the atom-independent 
continuum-continuum delay. We also show that the pres- 
ence of many resonances decaying into Ar + (3p 5 ) and 
Ar + (3s3p 6 ) in the studied energy region affects the re- 
sult of the measurements in ways that are sensitive to 
experimental parameters. In a future study we will ex- 
tend the treatment to include consistent ab initio calcu- 
lations of the single photon and laser-assisted photoelec- 
tron emission delays, in a similar way as done is Ref. |64j . 

However, it has been shown for the helium atom that 
dynamical screening effects caused by the polarization of 
the parent ion by the IR field can affect the delay of elec- 
tron emission |65j . The rich structure of the Ar + spec- 
trum at the relevant energies [52] suggests that such an 
effect could be significant in the present case. A straight- 
forward possibility to include those effects in the calcula- 



tion is to move to a time-dependent approach, including 
explicitly the shape of the IR pulse in order to account 
for the full many-body dynamics of the process, both on 
the relatively long time-scale of autoionization, and the 
sub-femtosecond scale set by the XUV-pulse. 

For further developments we note that the compact- 
ness of the MCHF and B-spline close-coupling approach 
presented here makes it well suited for time-dependent 
many-body calculations on systems in which exact cal- 
culations are untractable. Only a limited number of 
such calculations exist today, see for example [53] I5BT - 
I69j . all of them performed within the Time-Dependent 
R-matrix framework [SSI EH IZJJ- As a next step we 
will combine the present approach with a solution of the 
time-dependent Schrodinger equation along the lines of 
Ref. [?H [73] , where the calculation is made without com- 
plex scaling and the final physical result are extracted 
by projection of the finite-time wave packet on K-matrix 
scattering states. Absorbing boundaries are then used 
in order to avoid reflections on the walls of the computa- 
tional box in the case of long simulations. We foresee here 
that the common description of the structure part in the 
time-independent and time-dependent calculation will be 
of significant value for the tuning of the approximations 
regarding many-body effects that has to be introduced in 
the time-dependent approach. 
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Appendix A 

Here, we give the details of the calculations used to 
produce Fig. [S] The model is aimed at a qualitative 
description of the spectral features in the range of the 
experiments performed in Refs. [5] [S] . 

According to Ref. [ST], the most prominent features 
in the 3s photoionization partial cross-section and to- 
tal cross-section, come from the resonances of the type 
3p 4 nl[LSir]n'l' with nl — 3d, 4s, and Ap, n' < 7, and for 
the parent term LStt = 2 S e , 2 P°, 2 P e and 2 D e . The 
2 P e parent states are mostly responsible for resonances 
at 31 eV and below so we do not consider them. Us- 
ing the same Is, 2s, 2p, 3s and 3p orbitals as before, we 
perform Hartree-Fock level average calculations on the 
3p 4 nl states optimizing the 3d, 4s and 4p orbitals. Ex- 
tending the active set of spectroscopic orbitals by a 4d 
correlation orbital, we perform a MCHF calculation on 
the lowest state of the 2 S e symmetry of Ar + including all 
single and double (SD) excitations of the multi-reference 
set (MR) {3s3p 6 U3p 4 4sU3p 4 3dU3p 4 4d}. Configuration 
interaction calculations are then performed on the 2 P° 
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TABLE III. Comparison of the lowest computed levels of Ar^ 
to R-matrix calculations [5T] and experimental data |74) . 



State 


this work 


Ref. ED 


Exp. [74] 


3p b 2 P° 


16.481 


16.303 


15.819 


3s3p 6 2 S e 


29.24 


29.24 


29.24 


3p 4 ( 1 D)4s 2 D e 
3p i ( 3 P)3d 2 D e 
3p 4 ( 3 P)4p 2 P° 
3p 4 ( 1 S)4s 2 S e 
3p 4 CD)4p 2 P° 
3p i ( 1 D)3d 2 D e 
3p i ( 1 S)3d 2 D e 
3p 4 ( 1 D)3d 2 S e 


34.699 


34.999 


34.203 


35.418 


35.427 


34.462 


36.044 


36.514 


35.605 


37.394 


37.207 


36.504 


37.654 


38.134 


37.137 


38.647 


38.508 


37.148 


39.498 


39.296 


38.043 


39.595 


39.876 


38.585 



TABLE IV. Computed Fano profile parameters for the 
3s3p 6 np resonances. Positions are shifted upward by 0.15 eV. 



n 


AE 


r 


g(L) 


9(V) 




eV 


meV 






4 


26.524 


81.0 


-0.497 


-0.646 


5 


28.005 


25.0 


-0.435 


-0.489 


6 


28.516 


11.7 


-0.409 


-0.443 


7 


28.764 


6.4 


-0.396 


-0.423 


8 


28.903 


3.9 


-0.389 


-0.412 



and 2 D e states following the same MR-SD scheme with 
the references {3p 4 4sU3p 4 3dU3p 4 4d} and {3p 5 U3p 4 4j?}, 
respectively. Then, the CSF with a mixing coefficient 
over 0.07 are selected to build the target states. The 
lowest levels of Ar + obtained in this model are compared 
to experiment in Table |III| Most values agree with the 
experiment within 1 eV. 

The localized states Xf_ are chosen such that the close- 
coupling expansion Eq. ([5]) still retains the completeness 
of the B-spline set. Using this model, the photon en- 
ergy scale is shifted by 0.15 eV upward in order to match 
the experimental 3s photoionization. We choose a knot 
sequence as before for the inner region up to 1/Z ao, 
but with a linear step of 1 ao up to 51 ao and then an 
exponentially increasing step up to 150 a for accommo- 
dating orbitals with n < 7. The complex scaling starts 
at 40 ao with a complex scaling angle of 0.11 rad. As 
in Section \TV\ there is only a negligible difference be- 
tween the sum of the partial cross-sections and the total 
cross-section computed with Eq. (JsJ) . In Table IV we 
provide the results obtained for the 3s3p e np resonances 
line shapes. A good consistency is found with the val- 
ues given in Table [TTJ except for the q parameters. Fur- 
thermore, the partial cross-sections obtained within this 
model exhibit strong features in qualitative agreement 
with R-matrix calculations |61j . 
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